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ABSTRACT 

We obtain constraints on the variation of the fundamental constants from the fuh shape 
of the redshift-space correlation function of a sample of luminous galaxies drawn from 
the Data Release 9 of the Baryonic Oscillations Spectroscopic Survey. We combine this 
information with additional data from recent cosmic microwave background, baryon 
acoustic oscillations and Hq measurements. We focus on possible variations of the fine 
structure constant a and the electron mass me in the early universe, and study the 
degeneracies between these constants and other cosmological parameters, such as the 
dark energy equation of state parameter wde, the massive neutrinos fraction f^, the 
effective number of relativistic species iVeff, and the primordial helium abundance Yhb- 
In the case when only one of the fundamental constants is varied, our final bounds 
are a/ao = 0.9957!to;oo42 ^^^^ "ie/(me)o = 1.006tQ;gi[3. For the joint variation of both 
fundamental constants, our resuhs arc a/ao = 0.990lt^;^^^^ and me/(TOe)o = 1.028 ± 
0.019. The variations of a and me from their present values affects the bounds on other 
cosmological parameters. Although when me is allowed to vary our constraints on wde 
are consistent with a cosmological constant, when a is treated as a free parameter we 
find wde ~ —1.20 ± 0.13; more than 1 a away from its standard value. When fi, and 
a are allowed to vary simultaneously, we find < 0.043 (95% CL), implying a limit 
of < 0.46 eV (95% CL), while for me variation, we obtain < 0.086 (95% 

CL), which implies J] < 1.1 eV (95% CL). When N^s or Ihc are considered as 
free parameters, their simultaneous variation with a provides constraints close to their 
standard values (when the Hq prior is not included in the analysis), while when me is 
allowed to vary, their preferred values are significantly higher. In all cases, our results 
are consistent with no variations of a or rrie at the 1 or 2 cr level. 

Key words: cosmological parameters, large-scale structure of Universe, early Uni- 
verse 
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1 INTRODUCTION 

The overwhelming amount of cosmological observations ob- 
tained over the past few years has aUowed not only the pre- 
cise determination of the parameters of the standard cosmo- 
logical model but also has provided plenty of scope to test 
non-standard physics and cosmological assumptions such as 
the constancy of fundamental constants over cosmological 
timescales. 

The variation of fundamental constants is a predic- 
tion of theories attempting to unify the four interactions 
in nature, such as string derived field theories, r elated 
brane-world theories and Kal uza-Klein theories (see lUzanI 
l2003l : lGarci'a-Berro et al ] |2007l . and references therein) . Sub- 
stantial work have been devoted to constrain such varia- 
tions using cosm ol ogical observations (jRahmani et ahl 20121: 
Coc et all I2OIII: iLevshakov et al.l |2012| : iMenegoni et all 



20121 : iLandau fc Scoccolal I2OI0I ). Unifymg theories predic , 



the variation of all coupling constants, being all variations 
related in general to the rolling of a scalar field. In this pa- 
per we adopt a phenomenological approach and analyse the 
possible variation of the fine structure constant a and of the 
electron mass rrie between the recombination epoch and the 
present time, without assuming any theoretical model. 

The cosmic microwave background (CMB) is a pow- 
erful tool to study the early universe. The acoustic 
oscillations present in the CMB power spectrum are 
also imprinted through the bary ons on the large-scale 
structure (LSS) po wer spectrum (|Eisenstein fc Hul 1 19981 : 
iMeiksin et al.lll999l ). The correlation function ^(s) is the 
Fourier transform of the latter, and the oscillation struc- 
ture appears there as a single peak whose po sition is re- 
lated to the sound horizon at the drag redshift (|Matsubaral 
The o ngoing Baryonic Osc illation Spectroscopic Sur- 
vey (BOSS, iDawson et al.l I2OI2I') is a part of Sloa n Digi- 
tal Sky Survey-Ill (SDSS-IlT lEisenstein et~al]|201ll 'l and is 
aimed at obtaining redshifts for 1.5 x 10® massive galaxies 
out to z = 0.7 over an area of 10,000 deg^. BOSS is de- 
signed to measure the baryon acoustic oscillations (BAO) 
signal to probe the expansion history of the universe. This 
information places complementary constraints on the varia- 
tion of fundamental constants. A high redshift galaxy sam- 
ple from BOSS Data Release 9 (DR9), denoted CMASS, is 
constructed through a set of colour- magnitude cuts designed 
to select a ro ughly volume-limited s ample of massive, lumi- 
nous galaxies (|Eisenstein et al.]|201ll . Padmanabhan et al. in 
prep.). The clustering properties of the BOSS CMASS sam- 
ple have been analysed in detail in a recent series of papers 
llAnderson et al.ll2012l: iManera et al.ll2012l:lReid et aljlioiS : 



Ross et al.ll2012l: ISamushia et al.ll2012l : ISanchez et al.([2oT2 
Toieiro et al.ll2012l 'i 



The position of the peak in the correlation function of 
galaxies can place constraints on the variation of fundamen- 
tal constants. Moreover, the full shape of the correlation 
function provides additional information that can break de- 
generacies, since some parameters vary the full shape, while 
others affect only the position and height of the BAO peak. 
We use the full shape of the correlation function of BOSS- 
CM ASS galaxies presented in I Sanchez et al.l (|2012l ). in com- 
bination with CMB observations, to place constraints on the 
time variation of fundamental constants in the early uni- 
verse. We focus on possible variations in the fine structure 



constant, a, and the electron mass, me, at the recombina- 
tion epoch. Strictly speaking, the acoustic fluctuations in 
the baryons are fr ozen in at the drag e poch rather than 
at last scattering (|Hu fc Sugivamj Il996h . This dynamical 
decoupling of the baryons from the photons occurs never- 
theless near recombination. Therefore, we assume that the 
values of the fundamental constants are the same through- 
out this epoch, though they can differ from their current 
values. We analyze the degeneracies with the basic cosmo- 
logical parameters, as well as with others, such as the dark 
energy equation of state, the neutrino mass, the effective 
number of relativistic species, and the primordial helium 
abundance. 

Limits on the present rate of variation of a and 
fj. — rrie/mp (where mp is the proton mass) are pro- 



vided by atomic clocks ( Prestage et al. 



200 ll: iBize et al.l l2003l: 



Marion et al.l 



■19951 : Sortais et al 
20031 : iFischer et al 



2004 iPeik et all I2OO4I'). Data from the Oklo natura l fis 



sion reactor ( Damour fc Dvsonlll996l : Fuiii et al. |200G|) and 
half-lives of long lived P decayers (|01ive et al.l 12004 ) al- 
low to constrain the variation of fundamental constants 
at 2; ~ 1. Absorption systems in the spectra of high- 
redshift quasars put additional constraints at different red- 
shifts. The method is based on the measurement of the 
separation between spectral lines in doublets and multi- 
plets, whose dependence on the c onstants vary among dif- 
ferent species (see for example IWebb et al.l 1 19991. 1200 ll: 
Murphv et al] l2003l: lAgaforiova et al.l I2OIII : iKanekar et alj 



2012l : IWendt fc Molardl2012l '). Ahhough the limits imposed 



by CMB and LSS are less stringent than the previous ones, 
they are important because they refer to earlier times. 

The paper is organized as follows. In Section[2]we show 
how the correlation function depends on the values of a and 
rUe at the recombination epoch. We describe the datasets 
used to place constraints on the variation of the fundamental 
constants and the statistical method performed. The mod- 
elling of the correlation function is also summarized. Sec- 
tion |3] presents the results for different parameter spaces. 
Conclusions are outlined in Section |4l 



2 METHODOLOGY 

We performed a statistical analysis to constrain the variation 
of a and me at the recombination epoch, together with other 
cosmological parameters varied. In this Section we describe 
the datasets used to obtain our results. Then, we summa- 
rize the modelling of the correlation function, and present 
a brief explanation of why the correlation function is an ef- 
fective observable to constrain the variation of fundamental 
constants. Finally, we describe the statistical analysis em- 
ployed. 



2.1 Data 

In this paper, we use the full shape of the large- 
scale two-point correlation functio n ^(s) of the BOSS- 
CM ASS galaxy sample, computed in lSanchez et al.l (|2012h . 
This function was computed using the first spectro- 
scopic data release of BOSS (Data Release 9, DR9, 
ISDSS-III Collaboration: Ahn et al.ll2012h . The galaxy tar- 
get selection of BOSS is divided in two separate sam- 
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pies, named LOWZ and CMASS, covering different red- 
shifts ([Eisenstein et al ] |201ll : IPadmanabhan et alJIin prep. I : 
iDawson et all 120121 '). This selection is based on photomet- 
ric obs ervations done wit h the dedicated 2.5-m Sloan Tele- 
scope (|Gunn et al.l |2006| ), located at Apache Point Ob- 
servatory in New Mexico, usi ng a drift- scanning mosaic 
CCD camera ((Gunn et al.lll998l ) that produces ugriz images 
(|Fukugitaet al.lll996l ). The CMASS sample is constructed 
on the basis of gri colour cuts designed to select lumi- 
nous galaxies such that they constitute an approximately 
complete galaxy sa mple down to a limiting stellar mass 
iMarastonl lin prep.h . Spectra of the LO"WZ and CMASS 
samples are obtained using the double-armed BOSS spec- 
trographs, whi ch are significantly upgraded from those used 
by SDSS-I/II (|York et all l200of ). covering the wavelength 
range 3600 A to 10000 A with a resolving power of 1500 
to 2600 (|Smee et al.ll20l3 ). Spectroscopic redshifts are mea- 
sured usi ng the minimum-y^ template-fitting procedure de- 
scribed in lAihara et al.l (|201ll ). with te mplates and method s 
updated for BOSS dat a as described in lBolton et al.l (|2012t ). 
lAnderson et al.l (|2012l 'l present a detailed description of the 
construction of the catalogue for LSS studies based on this 
sample, and the calculation of the completeness of each sec- 
tor of the survey mask, i.e., the areas of the sky covered by 
a unique set of spectroscopic tiles. 

To constrain high-dimensional parameter spaces, it is 
necessary to combine the CMASS ^(s) with other datasets. 
In our analysis, we use the 'WMAP 7- year temperature 
and temperature-polarization power spectra (|Larson et al.l 
20 111') and the resu lts from the South Pole Telescope (SPT, 



Keisler et al.|[201ll ). The latter provide information on the 



structure of the acoustic peaks in the CMB po wer spectrum 
up to multipoles I ~ 3000. As discussed in iKeisler et al.l 
l|201ll) . in the multipole range covered by SPT (650 ^ ^ ^ 
3000), the CMB power spectrum contains a non- negligible 
contribution from secondary anisotropies, while for £ < 650 
the main contribution arises from primary a nisotropies. We 
follow the approach of iKeisler et al.l (|201ll ) to account for 
the secondary anisotropies, by including the contribution 
from the Sunyaev-Zel'dovich (SZ) effect and the emission 
from foreground galaxies (including both a clustered and a 
Poisson source contribution), using templates whose ampli- 
tudes are considered as nuisance parameters and marginal- 
ized over. The WMAP-SPT combination is referred to as 
our "CMB" dataset. 

We also use information from other clustering mea- 
surements in the form of constraints on the position of 
the baryon acous tic peak from independent analyses. We 
use the results of iBeutler et al] (|201ll ') which are based on 
measurements of the large-scal e correlation function of the 
6dF Galaxy Survey (6DFGS, Ijones et all 120091') and the 
2% di stance measurement ob tained by Padmanabhan et al.l 
(|2012h and IXu et all (|2012l ) from the application of an 
up dated version of t he rec onstruction technique proposed 
by lEisenstein et al.l (|2007l ) to the clus tering of g a laxies 
from the final SDSS -II LRG sample (|York et al'l l2000l : 
lEisenstein et al.l 1201 ll ). The results of these analyses are 
combined in the "BAO" dataset. 

Lastly, we consider a Gaussian prior on the Hubble pa- 
rameter based on the latest Hubble Space Teles cope (HST) 
obser vations of Ho = 73.8 ± 2.4 kms~^ Mpc~^ (|Riess et al.l 
H). 



The aforementioned datasets are used in different com- 
binations to check the consistency of the obtained bounds. 
Firstly, we use the CMB data alone, and then combine it 
with the CMASS correlation function. In the end, we com- 
bine the four datasets to obtain our final constraints. 

We do not consider supernovae (SNs) type la data be- 
cause the light curves of the SNs are obtained assuming 
that the fundamental constants have their present values at 
the observing redshift. However, since we are investigating 
a possible time evolution in the value of a and me, and 
the SNs are at consider ably high redshift (0.7 < z < 1.4 
for the high-z sample of lConlev et ahlboilh . we cannot ne- 
glect the possibility that the constants have a different value 
at those times. In fact, several studies aiming at measur- 
ing the value of a at high redshift using quasar absorption 
systems do not conclusively exclud e the variation of f unda- 
mental constants at those r edshifts ( Webb et al ] | 19991 . 12OIII : 
iMurphv et all |2003| . |2004| : iKing et all 120121 '). Therefore, to 
be conservative, we choose not to consider the supernovae 
datasets in our analysis. 



2.2 Model for the correlation function 

We follow ISanchez et all (|2012l ) and model the shape of the 
large-scale correlation function, (,{s), by applying the follow- 
ing parametrization: 



+ A 



MceUs)c«(.)], (1) 



where the symbol ® denotes a convolution, and the bias 
factor b, mode-coupling amplitude Amc, and the smoothing 
length fc* are considered as free parameters and marginal- 
ized over. Here is the derivative of the linear correlation 



function ^l, and ^l^'(s) is defined by 



C[^)(s) = s.V-^eL(.) = ^/ Pi 



(fc) ji (fcs)fc dfc, 



(2) 



with ji{y) denoting the spherical Bessel function of order 1. 

The parametrization o f equa tion [T] was first proposed 
by ICrocce fc Scoccimarrol (|2008l l ar id is based on renor- 
maliz e d perturbation theory (RPT, ICrocce fc Scoccimarrol 
bOOfih . I Sanchez et al.l (|2008l ) compared this model against 
the results of an ens emble of large volum e N-body simula- 
tions (L-BASICC-II. lAngulo et al.|[20o3 ). and showed that 
it provides an accurate description of the full shape of the 
correlation function, including also the effects of bias and 
redshift-space distortions. This parametrization has been 
applied to obtain constraints on cosmological parameters 
from clustering mea s urements from various galaxy sample s 
dSanchez et al.ll2009l: iBeutler et al.ll201ll : iBlake et al.]|201ll ') . 



As in ISanchez et all (|2012f ). we restrict the comparison 
of the model of equation [T] and the BOSS-CMASS corre- 
lation function to 40 < s < 200/i^^Mpc, and assume a 
Gaussian likelihood function. 



2.3 Effects on the full-shape of ^(s) 

During the recombination epoch, the ionization fraction is 
determined by the balance between photoionization and re- 
combination. The most important effects of changes in a 
and rUe during this epoch are due to their influence upon 
Thomson scattering cross section ar ~ /Smlc^ and 
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the binding energy of hydrogen Bi = ^a^nieC^ ■ The ioniza- 
tion history is more sensitive to a than to me because of 
the Bi dependence on these constants. The main result is 
the shift of the epoch of recombination to higher z as a or 
me increases, which corresponds to a smaller sound horizon. 
These effects are imprinted into the matter power spectrum 
through the transfer function. Hence, they also affect the 
galaxy correlation function. Consequently, if at recombina- 
tion a or me had a value higher than the present one, the 
peak in the correlation function would appear at smaller 
scales, since the position of the peak is related to the size of 
the sound horizon at the drag epoch. 

Fig. [1] shows the effect of a variation in a and in me, 
during the recombination epoch, that we should expect in 
the correlation function of CMASS galaxies, in a flat uni- 
verse with cosmological parameters (cjb, ti^dm, t, h, Us) = 
(0.0221, 0.1145, 0.696, 0.962). We show the prediction for 
different values of the fundamental constants: their present 
value, and a variation of ±5% with respect to their values to- 
day. Data points are the spherically averaged redshift-space 
two-point c orrelation function of the full CMASS sample 
presented in ISanchez et ai] (|2012l '). 

The changes in the correlation function due to varia- 
tions in a are larger than the changes in me for a given 
relative variation of their value. Both constants affect the 
position and height of the peak, leaving the rest of the 
curve unchanged. This effect can break degeneracies with 
other cosmological parameters that affect the full shape of 
the correlation function, such as the dark energy equation 
of state. Other parameters also affect the position of the 
peak. This leads to degeneracies between the fundamental 
constants and other cosmological parameters, such as fim, 
as we will see in Sections 13.11 and 13.21 The combination of 
CMASS data with different datasets helps to break those 
degeneracies. 



2.4 Statistical Analysis 

We perform our statistical analysis by exploring the pa- 
rameter spaces with Mont e Carlo Markov chain s generated 
with the COSMOMC code (iLewis fc Br idle 2002'). which uses 
the Boltzmann cod e CAMB (|Lewis et al..,200Q ) and recfast 
jSeager et all Il999l ) to compute the CMB power spectra. 
In order to be able to study more general models in which 
the dark energy component is different from the cosmolog- 
ical constant, we use a generalized version of CAMB which 
supports values of the dark energy equ ation of state be - 
yond the phantom divide, wde < —1 (|Fang et al.l |200§ ). 
We modified these codes to include t he variation in a an d 
rUe at recombination as described in iLandau et all (|2008h . 
Additional modifications from iKeisler et al.l ( 201ll ) are in- 
cluded to compute the likelihood of the SPT dataset. The 
dependence on the fundamental constants of the detailed 
ph ysics relevant in the recombination process is described 



Scoccola et al.l (|2008l ). Nevertheless, we emphasize that 



such description is done in terms of a modification to the 
effective 3-level atom model which is used in Recfast 
JWong et al. 200 §). Additional physical proc esses (see e.g., 
iRubifio-Martm et al.|[200i : iFendt et al.ll2009l . for a review) 
are effectively t reated using a correction function inside 
Recfast vl.5 (|Rubifio-Martm et al]|2010l '). This function 
has, to first order, a negligible dependence on the standard 



0.05 



0.04 



0.03 



0.02 



0.01 



BOSS-CMASS 
a/ag= 0.95 
alaQ= 1 .00 
oJaQ= 1 .05 



iJiiiiiiiiii^ 



50 



100 
s [h'"' Mpc] 



150 



200 



0.05 



0.04 



0.03 



0.02 



0.01 



BOSS-CMASS 

"'4'^oh= 0-95 
me/(mg)o=1.00 
mg/(mg ' " 



1.05 



50 



100 150 
s [h'"' Mpc] 



200 



Figure 1. Effects on the two-point correlation function of a 5% 
variation in a (upper panel), and in rrte (lower panel) at recom- 
bination time, with respect to their present values. Other param- 
eters are kept fixed 

. Data points are the measure ment of the BOSS-C MASS two- 
point correlation function from lSanchez et al. I HoHl). 

cosmological parameters , but also on non-stan dard parame- 
ters such as Vhc or A'efi (jShaw fc Chluballionl ). Although it 
has not been demonstrated explicitly, this correction func- 
tion is also expected to have a small dependence on a and 
me. Thus, we will use in this paper the correction function 
as it appears in Recfast vl.5. However, for future CMB 
experiments with higher sensitivities in the damping tail of 
the angular power spectrum, a more detailed and complete 
treatment of the recombi nation problem might be relevant 
(|Chluba fc Thomas|[20Tll ). 

We consider a spatially-flat cosmological model with 
adiabatic density fluctuations, described by the following 
parameters, which define the ACDM model: 

P = (cJb,i^dm, 6,T, y4s,71s) (3) 

where a;b ~ flhh^ is the baryon density, ojdm ~ ^dmh^ is the 
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dark matter density, both in units of the critical density; 
& gives the ratio of the comoving sound horizon at decou- 
pling to the angular diameter distance to the surface of last 
scattering; r is the reionization optical depth; and As and 
Us are the amplitude and spectral index of the primordial 
power spectrum of the scalar fluctuations, respectively, at 
the pivot wavenumber of fc = 0.05 Mpc~^. 

To this set of cosmological parameters, we add the value 
of a and/or me at the recombination epoch as additional 
variables. We introduce them as a/ao and me/(me)o, i.e., 
relative to their present value, denoted with a subscript. 

Moreover, to constrain possible deviations from the 
ACDM model, we study also the cases in which we allow 
for variations of the dark energy equation of state uide and 
the dark matter fraction in the form of massive neutrinos, 
= Qi^/ildm, focusing on the degeneracies with the fun- 
damental constants studied here. In addition, we analyze 
the case in which the effective number of relativistic species 
is different from its standard value of A^eff ~ 3.046, and 
the case in which the primordial helium fraction can have 
a value different than the standard one of Yuc = 0.24. Fi- 
nally, we present the constraints on the parameters of the 
ACDM-f iVeff and ACDM-fFHc models in the case of no vari- 
ation of the fundamental constants. 

We also present constraints on other quantities, derived 
from the above parameter set. These are the sum of the 
neutrino masses, given by 

^ m„ = 94.4 U eV, (4) 

the dark energy density Ode, the total matter density fim, 
the rms linear perturbation theory variance in spheres of 
radius 8/i~^Mpc erg, the age of the universe to, the redshift 
of reionization z^c, and the bubble factor h. 

The parameters b, Amc, and fc* of the model for the 
full shape of the correlation function, and the amplitudes 
Dfooo; ^3(K)0) ^3000 1 necessary to account for the secondary 
anisotropics in the SPT power spectrum, are treated as nui- 
sance parameters, and marginalized over when presenting 
our results. 

For most of the cases, we only use the information from 
the CMASS correlation function in combination with our 
CMB data. However, in Sections 13.11 and 13.21 we also ran 
chains using the CMASS information alone, to study the 
degeneracies between a and rrie and other cosmological pa- 
rameters obtained from this dataset. In these cases we im- 
pose Gaussian priors on ujh and Us, obtained from CMB-only 
results (see first column of Table [1] and Table [2]) . 



3 RESULTS 

We present the constraints obtained for the fundamental 
constants and cosmological parameters in each of the cases 
studied. Errorbars will indicate 68% confidence level (CL) 
unless otherwise stated. In Section lSTll we study the variation 
of the fine structure constant, and do the statistical analy- 
sis varying also the cosmological parameters of the ACDM 
model. In Section [3]2] we study the variation of the electron 
mass together with the cosmological parameters. Section [53] 
investigates the joint variation of a and rUe . In Section 13.41 
we analyze the constraints in the case where one of the con- 
stants (a or rrie) and the dark energy equation of state wde 



can take values that differ from the standard ones. In Sec- 
tion 13.51 we study the case where we vary one of the fun- 
damental constants and the massive neutrinos fraction /j^. 
Section 13.61 focusses on the constraints obtained when the 
fundamental constants are varied together with the effec- 
tive number of relativistic species, A'efr, while Section 13.71 
presents the constraints when the fundamental constants are 
varied together with the primordial helium fraction Yue- In 
Appendix [XI we give the constraints on the ACDM model 
for our four datasets, to facilitate the comparison to the re- 
sults presented in this paper. 

3.1 Variation of a 

In this Section, we extend the ACDM model to include pos- 
sible variations in the fine structure constant during recom- 
bination. We present our constraints on a/ao and the cos- 
mological parameters. 

The upper panel of Fig. [2] shows the two-dimensional 
marginalized constraints in the a - Sim plane obtained from 
the CMB and CMASS datasets in isolation. As specified in 
Section [2j when using the information of the CMASS corre- 
lation function alone we impose Gaussian priors on toh and 
ris consistent with our CMB-only results. The constraints 
obtained from these datasets exhibit strong degeneracies. 
As these degeneracies constrain different combinations of 
a and f2m, the combination of the two datasets provides 
tighter constraints on both parameters simultaneously. As 
can seen in the lower panel of Fig. [2j the remaining de- 
generacy between these parameters is alleviated when more 
datasets are added to the analysis. The correlation factors 
are -0.62, -0.38, and -0.13 for CMB, CMB+CMASS, and 
the full dataset combination (CMB+CMASS+BAO+i^o), 
respectively. 

In Fig.|3]we show the constraints in the a - Ho plane. Ho 
is better constrained when additional datasets are included 
in the analysis. The value of the HST prior on Ho lies almost 
2 a from the value obtained for Ho from the CMB-I-CMASS 
dataset, being the latter considerably smaller. Hence, when 
the HST prior on Ho is taken into account, the obtained 
value for the Hubble parameter is increased. Due to the de- 
generacy with a, the contours are shifted in the parameter 
space when the Ho prior is included, and the value of the fine 
structure constant is increased, being closer to its present 
value. 

In Table [1] we present the constraints on the variation 
of a and the cosmological parameters. The constraint on 
a from CMB data alone is a/ao = 0.9914 ± 0.0055. The 
CMB-only constraints on Qrn are considerably degraded by 
the inclusion of a as a free parameter, with fim ~ 0.3031q o36 
(compare to the ACDM value of Qm in Table IXT1 of the Ap- 
pendix). When we add the information encoded in the full 
shape of the CMASS ^(s), we find a/ao = 0.9917 ± 0.0046 

and fim = o.295iE;;;;i^. In the case that all datasets are used 

in the analysis, the constraints are a = 0.99571q oo42 and 
Qm = 0.283 ± 0.010, which is the same precision as that ob- 
tained in the ACDM model. The inclusion of the additional 
datasets produces only a mild improvement of the bounds 
on a with respect to those obtained using CMB information 
alone. This result arises because our CMB dataset covers the 
high multipole range, with £ ~ 3000, which imposes strong 
constraints on this parameter. The bounds obtained for the 
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Figure 2. The marginalized posterior distribution in tlie a - Qrn 
plane for the ACDM parameter set extended to include the vari- 
ation of a. The dot-dashed lines show the 68% and 95% contours 
obtained using CMB information alone. The dotted contours, in 
the upper panel, show the results from CMASS correlation func- 
tion. The dashed lines correspond to the results obtained from 
the combination of CMB data plus the shape of the CMASS 
^(s). The solid lines indicate the results obtained from the full 
datasct combination (CMB-t-CMASS-t-BAO-t-i^o)- The vertical 
dotted line corresponds to the ACDM model, with q/qq = 1 



cosmological parameters from the full dataset are consistent 
within 1 a with their values in the ACDM model (see Ta- 
ble ET^ 



Menegoni et~al] (|2012l ') present constraints on the vari- 



ation of Q using CMB data including data from SPT and 
the A tacama Cosmology Telescope (ACT, iDunklev et al.l 
I2OIII ). both probing the damping regime of the CMB fluctu- 
ations. By combining this information with the galaxy power 
spec trum from the S DSS-DR7 luminous red galaxy sam- 
ple (|Reid et al.lbOld 'l and the HST prior on Ho, they find 
a/ao — 0.984 ±0.005. The precision on this constraint on a 



Figure 3. The marginalized posterior distribution in the a - 
Ho plane for the ACDM parameter set extended to include the 
variation of a. The dot-dashed lines show the 68% and 95% 
contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination 
of CMB data plus the shape of the CMASS $(s). The solid lines 
indicate the results obtained from the full datasct combination 
(CMB+CMASS+BAG+i^o)- 



is provided mostly by the CMB data, explaining why our re- 
sults with CMB data alone have almost the same precision 
as theirs. Nevertheless, the CMASS ^(s) further improves 
the precision of the bound. It is important to emphasize, 
however, that our results prefer values for a that are closer 
to its present value than theirs. Our results are consistent 
with no variation of a at 2 a (for the full dataset, they are 
just slightly inconsistent at 1 a). 



3.2 Variation of rUe 

Now we turn to the case in which we extend the ACDM 
model to include possible variations in the electron mass 
during recombination. 

The contours in Fig. |3] show the two-dimensional 
marginalized constraints in the me - fim plane. The upper 
panel shows that the constraints obtained from the CMB 
and CMASS datasets in isolation exhibit strong degenera- 
cies in different directions in the parameter space. The CMB- 
only constraints exhibit a strong degeneracy between rUe and 
f2m, causing Q^n to be poorly constrained by the CMB data 
alone, with = OM^o ll- The inclusion of the CMASS 
correlation function improves the constraints on by a 
factor larger than five to obtain Qui = 0.295 ± 0.021. This 
bound is further reduced by a factor of two when all datasets 
are included, in which case we find Qui = 0.280 ±0.010, hav- 
ing the same precision than in the ACDM model. 

In Fig. [S] we show the two-dimensional marginalized 
constraints in the - Ho plane. When rrie is allowed to 
vary, CMB information alone is insufficient to place any re- 
liable constraint on Ho- When the CMASS dataset is added, 
the constraint improves noticeably. This dataset breaks the 
degeneracy between Ho and m^. When all the datasets are 
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Table 1 . The marginalized 68% allowed regions on the cosmolog- 
ical parameters of the ACDM model, adding the variation of the 
fine structure constant, a, obtained using different combinations 
of the datasets. 
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+ BAO + f^o 
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considered, the value of Ho is increased due to the HST prior, 
and the value of me is shifted towards its present value. 

In Table[2]we present the constraints on the variation of 
me and the cosmological parameters. The constraint on me 
from CMB data alone is 0.989t55 g^^. When we add the in- 
formation encoded in the full shape of the CMASS £,{s), the 
bound is 0.98llo o2i- the case that all datasets are used 
in the analysis, the constraint is 1.006to'oi3- The precision 
in the bound is highly improved when we add the informa- 
tion of the CMASS ^(s) to the analysis. Our final results are 
completely consistent with no variation of nie within 1 a. 

The bounds for the cosmological parameters are consis- 
tent within 1 a with those from the ACDM model. The mean 
values of as and h for the CMB+CMASS dataset are some- 
what smaller than in the ACDM model, but still consistent 
within 1 (7. 



3.3 Joint variation of a and me 

In this Section, we extend the ACDM model to study the 
joint variation of a and me. 

The contours in Fig. [6] show the two-dimensional 
marginalized constraints in the me - ce plane. When the Hq 
prior is used, the mean value of a is marginally decreased, 
while the mean value of me is increased. The inclusion of 
additional datasets reduces the allowed region in the param- 
eter space, while increasing the correlation between the fun- 
damental constants; the correlation factor is —0.23, —0.55, 
and -0.68 for CMB, CMB+CMASS, and the fuU dataset 
combination (CMB+CMASS-|-BAO-|-/fo), respectively. 

In Table [3] we present the constraints obtained for the 
fundamental constants and the cosmological parameters. 
When adding different datasets, the precision in the deter- 
mination of a remains the same, although the mean value 
is slightly lower when the full dataset is used. For the CMB 
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Figure 4. The marginalized posterior distribution in the rUe 
- Qm plane for the ACDM parameter set extended to include 
the variation of me. The dot-dashed lines show the 68% and 95% 
contours obtained using CMB information alone. The dotted con- 
tours, in the upper panel, show the results from CMASS correla- 
tion function. 

The dashed lines correspond to the results obtained from the com- 
bination of CMB data plus the shape of the CMASS The 
solid lines indicate the results obtained from the full dataset com- 
bination (CMB+CMASS+BAO-f-i^o)- The vertical dotted line 
corresponds to the ACDM model, with r?ie/(rrae)o = 1- 



and CMB-I-CMASS datasets, the mean value of a is almost 
the same, regardless of me being fixed to its present value 
or allowed to vary. For the full dataset, the value of a is de- 
creased by 1 (J when rUe is also allowed to vary (see Fig. [7]) . 

The constraints on me are significantly improved when 
adding extra datasets. The bound on me from CMB-only 
data is improved by almost a factor of three when the 
CMASS dataset is added to the analysis, and is further 
improved when all the datasets are considered. The mean 
value of me is increased when a is also allowed to vary, com- 
pared to the case in which a is fixed to its present value. For 
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Figure 5. The marginalized posterior distribution in the me - 
Ho plane for the ACDM parameter set extended to include the 
variation of me/(me)o- The dot-dashed lines show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS ^(s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO+/i'o)- 



Figure 6. The marginalized posterior distribution in the a - rUe 
plane for the ACDM parameter set extended to include the joint 
variation of a and rrie. The dot-dashed lines show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS £,[s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO+i^o)- 



Table 2. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, adding the variation of 
the electron mass, me, obtained using different combinations of 
the datasets. 





CMB 


CMB -1- CMASS 


CMB -1- CMASS 
+ BAG -f- Ho 


me/{me)o 


n QSQ+0 067 
0.989_o 069 




1 006+°°i'* 

i.UUD_o 013 


1000 


-, p,oo+0.049 
l.UJO_g 051 


J^-UZ/_o 015 


1 n4/lS+0-0097 
i.0448_o.oo94 


lOOiJb 


2.2olg:r7 


r, 1 77+0.054 
^■J^' '-0.054 


OOS+0.042 
z.zzo_o 042 


lOOl^dm 




lU.90_o g4 


^^•°'-0.54 


r 


0.0850l«:»?? 


0.0850^0.0069 


n nsi Q+0.0059 

0.08iJ_o.0066 


ns 


0.965lHi^ 


9646+° °"° 
U.yD1D_o 0099 


9620+° °°''* 

U.9DZU_o 0098 


In(lOiOAs) 


3.080tH^l 


p,7o+0.031 
o.u/o_o 030 


q p,oQ+0.029 
o.uo»_o 029 




69+° " 


U. 'U0_o 021 


0.720l°:°i° 






U-31_o.ii 


295+0 021 
U.^90_0 021 


0.280l°:°l° 




0.7991°:°?^ 


0.794t°:°l 


0.839t°:°l? 


to/Gyr 




14.20t°:f7 


13.64l°:l 




lo.slli 




lo.stiJ 


h 




0.668l°:°g 


704+°-0" 



the CMBH-CMASS and the full datasets, the increment is 
more than 1 a (see Fig. [H] for the full dataset case). For 
the full dataset combination, the mean value of me takes its 
largest value. Our final bounds, using the full dataset, are 
a/ao = 0.990lt|;;!^E;^4 and me/(me)o = im9,t°oMi9- Both 
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Figure 7. The marginalized, one-dimensional likelihood distri- 
bution for the fine structure constant, when only a is allowed to 
vary, and in the case of the joint variation of a and rrie- 

limits are consistent with no variation of the fundamental 
constants at the 2 a level. 

Most of the cosmological parameters are consistent 
within 1 a with their ACDM model values for each of the 
datasets. The value of ris is larger than in the ACDM model, 
for all of the datasets. For the CMB-I-CMASS and the full 
datasets, ris is larger by slightly more than 1 a. The value 
of (78 is slightly larger than its value in the ACDM model, 
but consistent within 1 o. 
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Table 3. The marginalized 68% allowed regions on the cosmolog- 
ical parameters of the ACDM model, adding the variation of the 
fine structure constant a and of the electron mass me , obtained 
using different combinations of the datasets. 
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Figure 8. The marginalized, one-dimensional likelihood distri- 
bution for the electron mass, when only me is allowed to vary, 
and in the case of the joint variation of a and me. 

3.4 Variation of fundamental constants and woe 

Until now, we have assumed that the dark energy compo- 
nent corresponds to a cosmological constant, with a fixed 
equation of state specified by uidb = —1. In this Section 
we explore the constraints on the value of woe, assumed to 
be redshift-independent, in the context of the variation of 
fundamental constants. 

In this study, the dynamical dark energy models are al- 
lowed to cross the so-called phantom divide, woe = —1, to 
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Figure 9. The marginalized posterior distribution in the a - 
Woe plane for the ACDM parameter set extended to include the 
variation of a and the redshift-independent value of wide as an 
additional parameter. The dot-dashed linos show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS ^(s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO-|-//o). The vortical and horizontal dotted 
lines correspond to the ACDM model, with a/ao = 1 a-nd '^DE = 
-1. 



explore models with wde < —1. In the framework of general 
relativity, a single scalar field cannot cro ss this threshold , 
since it becomes gravitationally unstable (|Feng et al.ll2005l : 
|Vikman.2005;.Hu..2005. ; ,Xia ct al...2008i ). Thus, more degrees 
of freedom are required, which are difficult to implement 
in general dark energy studies. We f ollow the pa r ametr ized 
post-Friedmann (PPF) approach of iFang et al] (|2008l ') , as 
implemented in CAMB, which provides a simple solution to 
these problems for models in which the dark energy compo- 
nent is smooth compared to the dark matter. 

First, we perform a statistical analysis varying a and 
Woe, together with the rest of the cosmological param- 
eters. Fig. |9] presents the constraints on the a - woe 
plane. Again, the inclusion of additional datasets reduces 
the allowed region in the parameter space, while increas- 
ing the anti-correlation between the fundamental con- 
stants. The correlation factor is —0.16, 0.37, and 0.54 
for CMB, CMB-fCMASS, and the full dataset combina- 
tion (CMB-^CMASS+BAO-h//o), respectively. For the fuU 
dataset, the value woe ^ — 1 is excluded at more than 1 a. 

Table|3]shows the constraints obtained for a, woe, and 
the cosmological parameters. The inclusion of more datasets 
slightly improves the constraint on a and does not appre- 
ciably affect the mean value. In the case of woe, a large 
improvement is observed, and the mean value of woe shifts 
towards lower values. Our final bounds, using the complete 
dataset, are a/ao = 0.9915±0.0048 and woe = -1.20±0.13. 
Our results are consistent with no variation of a at the 2 a 
level, while u;de is compatible with a cosmological constant 
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Table 4. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, adding the variation of 
the fine structure constant a, and the dark energy equation of 
state uijjEi obtained using different combinations of the datasets. 
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also within 2 a, although at 1 a, both quantities differ from 
their standard values. 

There is a slight tension at the 1 a level in the values of 
some of the cosmological parameters with respect to their 
values in the ACDM model. At 2 a, however, the results are 
consistent. Indeed, when compared to the ACDM model, we 
find that the value of O is lower and the age of the universe 
is higher for all of the datasets. The value of Us is higher for 
the CMB and the CMB+CMASS datasets, and as is higher 
for the full dataset. 

We then perform a statistical analysis varying and 
wuE- Fig. 1101 presents the resulting constraints on the rrie - 
wuE plane. CMB data alone cannot place strong constraints 
on the value of me and wue- The inclusion of the CMASS 
correlation function restricts the allowed region in the pa- 
rameter space, reducing the allowed range of values of me 
by a factor of two. Using the full dataset, our bounds are 
me/(me)o = 0.996 ± 0.029 and wde = -1.12 ± 0.23. In 
Table [5] we also present the constraints on the other cosmo- 
logical parameters. The standard case of me — (me)o and 
Wde = — 1 is completely consistent with all of the datasets. 
The values of the cosmological parameters are consistent 
within 1 a with their values in the ACDM model. 



3.5 Variation of fundamental constants and f,^ 

In the standard ACDM scenario, the dark matter component 
is given entirely by cold dark matter. However, neutrino os- 
cillations found in recent experiments imply that they have 
a non-zero mass that contributes to the total energy bud- 
get of the universe. A variation in the neutrino mass can 
alter the redshift of matter-radiation equality, thus mod- 
ifying the CMB power spectrum. Furthermore, until they 
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Figure 10. The marginalized posterior distribution in the me 
- «>DE plane for the ACDM parameter set extended to include 
the variation of me and the redshift-independent value of uide sts 
an additional parameter. The dot-dashed lines show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS ^(s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO+i?o). The vertical and horizontal dotted 
lines correspond to the ACDM model, with me/(me)o = 1 and 
Wue = — 1. 



Table 5. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, adding the variation of 
the electron mass me, and the dark energy equation of state iodEi 
obtained using different combinations of the datasets. 
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Figure 11. The marginalized posterior distribution in the a - fi, 
plane for the ACDM parameter set extended to include the vari- 
ation of a and a non- negligible fraction of massive neutrinos. The 
dot-dashed lines show the 68% and 95% contours obtained using 
CMB information alone. The dashed lines correspond to the re- 
sults obtained from the combination of CMB data plus the shape 
of the CMASS £,{s). The solid lines indicate the results obtained 
from the full dataset combination (CMB-|-CMASS-(-BAO-|-i?o)- 



become non-relativistic, neutrinos free-stream out of den- 
sity perturbations, suppressing the growth of structures on 
scales smaller than the horizon at that time, which depends 
on their mass, thus affecting the shape of the matter power 
spectrum and the correlation function. In this Section we 
explore the constraints on the neutrino fraction, when 
one of the fundamental constants (o? or me) is allowed to 
vary. We assume three neutrino species of equal mass. 

Fig. [11] shows the contours in the a - plane. The 
constraints on a are only slightly poorer than in the = 
case. When additional datasets are incorporated to the fit, 
the bounds on a tighten. 

In Table [6] we present the constraints obtained for a, 
fu , and the remaining cosmological parameters. In the case 
of only CMB data, we find < 0.098 at 95% CL, and 
a/ao — 0.9933j;Q oo5g . When we also include the informa- 
tion from the CMASS correlation function, this limit is re- 
duced to /i, < 0.049 at 95% CL, and the constraint on a is 
a/ao = 0.9940lo:oo49- The constraints for the full data set 
are < 0.043 at 95% CL and a/ao = 0.9978to:oo45- We 
find no degeneracy between a and When all datasets are 
included in the analysis, the value of a is increased (due to 
the larger value of h), our results being consistent with no 
variation of a within 1 a. 

Regarding the constraints on the primary cosmologi- 
cal parameters, we find that they are consistent with the 
ACDM values at 1 cr in most of the cases. For the CMB and 
CMB-I-CMASS datasets, some of the derived parameters dif- 
fer by more than 1 a from their values in the ACDM model, 
but they are consistent within 1 a when we consider the full 
dataset, with exception of as, which is lower by more than 
1 a also for the full dataset. For the CMB dataset, h differs 
by more than 2 a from the ACDM value, by 1 cr for the 



Figure 12. The marginalized posterior distribution in the nie - 
/,y plane for the ACDM parameter set extended to include the 
variation of me and a non-negligible fraction of massive neu- 
trinos. The dot-dashed lines show the 68% and 95% contours 
obtained using CMB information alone. The dashed lines cor- 
respond to the results obtained from the combination of CMB 
data plus the shape of the CMASS ^(s). The solid lines in- 
dicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO-l-//o)- 



CMB-I-CMASS dataset, and is fully consistent for the full 
dataset. The constraints at 95% CL on the sum of the neu- 
trinos masses are < 1.2 eV (CMB), '^m,^ < 0.53 eV 
(CMB-f CMASS), andJ2m^ < 0.46 eV (full dataset). 

Fig. [12] shows the contours for the analysis in the rUe - 
plane. The precision in the constraints on rUe are similar 
as those in the fv = case. When the CMASS correlation 
function is added to the CMB dataset, the allowed region 
in the parameter space is reduced. Due to the degeneracy 
between Ho and rrie, and between the latter and fi, (corre- 
lation factor equal to 0.83 and 0.72, respectively), when the 
full dataset is considered, the contours are shifted towards 
a region with higher values of rrie and 

In Table [7] we present the constraints obtained for me, 
/y, and the cosmological parameters. In the case of only 
CMB data, we find < 0.11 at 95% CL, and me/(me)o = 
l.Olllg 066- When we also include the information from the 
CMASS correlation function, this limit is reduced to < 
0.074 at 95% CL, and the constraint on rrie is me/(me)o = 
1.0091 

o!o28- The constraints for the full data set are < 
0.086 at 95% CL and me/(me)o = 1-035 ± 0.021. When 
all datasets are included in the analysis, the value of rrie is 
increased (due to the larger value of h), and our results are 
consistent with no variation of rrie at 2 cr. 

With regards to the constraints on the cosmological pa- 
rameters, we find that they are consistent with the ACDM 
values at 1 a in most of the cases, although ujdm differs from 
the ACDM value by 1 cr for the full dataset. There is also 
tension for erg at the 1 a level for all of the datasets, in the 
sense that our constraints are lower than the ACDM val- 
ues. For the full dataset, the age of the universe is lower 
and h is higher by more than 1 a than their corresponding 



© 0000 RAS, MNRAS 000, 000-000 



12 



C.G. Scoccola et al. 



Table 6. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model, adding the variation of the fine 
structure constant a and fi, as a free parameter, obtained using different combinations of the datasets. 
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Table 7. The marginalized 68% allowed regions on the cosmological parameters of the ACDM model, adding the variation of the 
electron mass rrte and fy as a free parameter, obtained using different combinations of the datasets. 
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values in the ACDM model. The constraints at 95% CL on 3.6 Variation of fundamental constants and A^off 

the sum of the neutrinos masses are < 1.4 eV (CMB), 

^rriiy < 0.91 eV (CMB+CMASS), and '^m^ < 1.1 eV (full The effective number of relativistic species, A'^oft, has been 
dataset). reported as higher than the expected standard va lue of 

A^cff = 3.046 (jPunklev et al.ll201l| - IXeisler et al.ll201ll ). This 
may indicate either additional relativistic species, or evi- 
dence for non-standard decoupling. Here we study the con- 
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Figure 13. The marginalized posterior distribution in the a 
- N^ff plane for the ACDM parameter set extended to include 
the variation of a and a variable effective number of relativistic 
degrees of freedom Ngg. The dot-dashed lines show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS ^(s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO+i^o). 



straints from the CMASS correlation function on the vari- 
ation of a and rrie when A'cff can differ from its standard 
value, providing details on the constraints set by the differ- 
ent datasets considered in this paper. 

In the case of no vari ation of the fundamental con- 
stants, iKeisIer et all (|2oi3) present a bound of A'cft = 
3.85 ± 0.62 for the CMB dataset (WMAP+SPT). We con- 
sider the ACDM-|-A^off model and obtain the constraints 
A^cff = 3.75 ± 0.58 for CMB-hCMASS and iVefi = 3.86±[;:;i^^ 
for the full dataset. 

The contours in Fig. [TH] show the two-dimensional 
marginalized constraints in the NcB - a plane. There is a 
strong degeneracy between a and N^a in the CMB dataset. 
When we add the CMASS dataset, this degeneracy is par- 
tially alleviated, and it disappears for the full dataset. The 
inclusion of the HST Hq prior shifts the contours towards 
higher values of A'cEF and a, making the latter consistent 
with its present value at 1 a. 

In Table[8]we present the constraints on A^ofi, ce, and the 
rest of the cosmological parameters for all of our datasets. 
When a is allowed to vary, N^s = 3.37toly for the CMB 
dataset, which is a value more consistent with the stan- 
dard one than when a is fixed to its present value (as 
said before, for WM AP-I-SPT data, the value presented by 
iKeisler et al]|201ll . is A^^^ff = 3.85 ± 0.62). The constraint 
on a is a/ao = 0.9935 ± 0.0093. By adding the CMASS 
correlation function, we tighten the bounds on NcU and q, 
which become N^s = 3.2lj;^;®? and a/ao = 0.9926 ± 0.0058. 
For CMB-hCMASS+BAO, A^ff = 3.19 ± 0.64 and a/ao = 
0.9926lg Q055 . When all datasets are considered in the anal- 
ysis, the constraints are Nch = 3.83 ± 0.40 and a/ao = 
0.9967 ±0.0042. The Ho prior shifts Aefi towards higher val- 



Table 8. The marginalized 68% allowed regions on the cosmolog- 
ical parameters of the ACDM model, adding the variation of the 
fine structure constant a and A^eff a-s a free parameter, obtained 
using different combinations of the datasets. 
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ues, which lie almost 2 a from the standard value of 3.046. 
O n the other hand , a is consistent with no variation at 1 
a. iMenegoni et all (|2012l ) studied the variation of a when 
Aeff can have values different from the standard one, al- 
lowing also for variations in the primordial helium abun- 
dance, Yuc- Their constraints are a/ao = 0.990 ± 0.006 and 
Aeff = 4.101°:^^ The difference in the datasets is that they 
use the DR7 LRG power spectrum and add the ACT data 
to their analysis, presenting their results only for the full 
dataset. 

There is tension at the 1 a level in some of the cos- 
mological parameters with respect to their ACDM model 
values. For the full dataset, G and to are more than 1 a 
lower, and Ua, erg, h and cjdm are larger by more than 1 a. 
Q is also more than 1 a lower than its concordance value 
for the CMB-I-CMASS dataset. In the rest of the cases, the 
parameters are all consistent within 1 a with their ACDM 
values. 

We now analyze the results for the case in which nie 
and Aefi are allowed to vary. In Fig. [14] we show the two- 
dimensional marginalized constraints in the Nes - me plane. 
There is no degeneracy between rUe and NbB for the CMB 
and CMB-I-CMASS datasets. When the fuU dataset is con- 
sidered, a mild correlation between both quantities appears. 

In Table [9] we present the constraints on Nch, rUe, 
and the rest of the cosmological parameters for all of our 
datasets. For the CMB dataset, N^s = 3.86 ± 0.60, which is 
more than 1 a higher than the standard value. The con- 
straint on nie is me/(me)o ~ 0.988 ± 0.065. By adding 
the CMASS correlation function, we significantly tighten 
the bounds on me but only slightly improve the limits on 
Acff, which become Acff = 3.80 ± 0.57 and 171^/(711^)0 = 
0.981 ±0.020. For CMB±CMASS±BAO, N^s = 3.78 ±0.58 
and me/(me)o ~ 0.981 ± 0.018. When all datasets are con- 
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Figure 14. The marginalized posterior distribution in tlic 
- NcS plane for the ACDM parameter set extended to inelude 
the variation of and a variable effeetive number of relativistie 
degrees of freedom A^efF- The dot-dashed lines show the 68% and 
95% contours obtained using CMB information alone. The dashed 
lines correspond to the results obtained from the combination of 
CMB data plus the shape of the CMASS ^(s). The solid lines 
indicate the results obtained from the full dataset combination 
(CMB+CMASS+BAO+Ho). 



Table 9. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, adding the variation of 
the electron mass rrie and A^^ff as a free parameter, obtained using 
different combinations of the datasets. 
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sidered in the analysis, the constraints are iVoff = 4.04l°;l^ 
and me/{me)o = O.gSStU JJl^. The i?o prior shifts N^s to- 
wards even higher values, which are more than 2 a above the 
standard value of 3.046. On the other hand, rrie is consistent 
with no variation at the 1 a level for all of the datasets. 

Regarding the values of the remaining cosmological pa- 
rameters, they are all consistent within 1 a with their ACDM 
values for the CMB dataset. For the CMB-I-CMASS dataset, 
ris is higher than its standard value by more than 1 a. For 



the full dataset, ojdn 



erg, and h are higher than their 



standard values in more than 1 a, while to is lower in more 
than 1 o. 



3.7 Variation of fundamental constants and Yu^ 

Light nuclei begin to form in a process known as big bang nu- 
cleos y nthesis (BBN.lAlpher et allll948l : ISchramm fc Turneij 



1 19981 : ISteigmanI |2007| ). when the universe cools to T ~ 
0.1 Me V. We denote the primordial abundance (mass frac- 
tion) of *He as Yhg, which is a function of the baryon density 
and the expansion rate during BBN. The value of Vho can 
be estimated by the effect of helium on the CMB damping 
tail. Helium combines earlier than hydrogen, and thus more 
helium (at fixed baryon density) leads to fewer free electrons 
during hydrogen recombination. This, in turn, leads to larger 
diffusion lengths for photons and less power in the CMB 
damping tail. We study the constraints from the CMASS 
correlation function on the variation of a and me when 
can differ from its standard value of "KHe = 0.24. 

In the ca s e of n o variation of the fundamental constants, 
iKeisler et all (|201ll ') present a bound of Vho = 0.296 ±0.030 
for the CMB dataset (WMAP-I-SPT). We consider the 
ACDM-|-Yhc model and obtain the constraints = 0.297± 
0.030 for CMB+CMASS and Fho = 0.30lto;o29 for the fuU 
dataset. 

The contours in Fig. [T5] show the two-dimensional 
marginalized constraints in the Yhg - a plane. There is a 
strong degeneracy between a and Yhb in all of the datasets. 
When a is allowed to vary, CMB information alone is in- 
sufficient to place any reliable constraint on Ihc- When the 
CMASS dataset is added, the constraint improves notice- 
ably. The inclusion of the HST Ho prior improves further 
the constraints, and shifts the contours towards higher val- 
ues of Vhc and a, making the former measurement inconsis- 
tent with its standard value at 1 ct. 

In Table [10] we present the constraints on Ihc, ct, 
and the rest of the cosmological parameters for all of our 
datasets. When a is allowed to vary, Yhb = 0.26 ± 0.15 
for the CMB dataset, which is a value consistent with the 
standard one, contrary to the case when a is fixed to its 
pre sent value (for WM AP-I-SPT data, the value presented 
bv iKeisler et al.ll201ll . is Fhc = 0.296 ± 0.030). The con- 
straint on a is a/ao = 0.994l|5;225 • By adding the CMASS 
correlation function, we tighten the bounds on Ihc and a 
to Fhc = 0.249 ± 0.064 and a/ao = 0.99201 g'.UEl^a . When 
all datasets are considered in the analysis, the constraints 
are Y^e = 0.314 ± 0.043 and a/ao = 1.0023li;;i;51^i . The Ho 
prior shifts Yhc towards higher values, which are almost 2 a 
above the standard value of 0.24. On the other hand, a is 
consistent with no variation at 1 a. iMenegoni et all (|2012l ) 
studied the joint variation of a, NcS and Yhc; their con- 
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Figure 15. The marginalized posterior distribution in the a 
- ^Ho plane for the ACDM parameter set extended to include 
the variation of a and a variable primordial helium abundance 
Yfjo- The dot-dashed lines show the 68% and 95% contours ob- 
tained using CMB information alone. The dashed lines corre- 
spond to the results obtained from the combination of CMB 
data plus the shape of the CMASS ^{s). The solid lines in- 
dicate the results obtained from the full dataset combination 
(CMB-l-CMASS-l-BAO+i^o). 



straint is Yhg = 0.215 ± 0.096 but with A^efi higher by more 
than 3 a than its standard value. 

There is a slight tension at the 1 a level for the pa- 
rameter Us, which is higher than its ACDM value for the 
CMB-(-CMASS and the full datasets. In the rest of the 
cases, the parameters are all consistent within 1 a with their 
ACDM values. 

We now analyze the results for the case in which rUe and 
Yhc are allowed to vary. Fig. [16] shows the two-dimensional 
marginalized constraints in the Yhc - rUe plane. There is no 
degeneracy between rUe and Yhg for any of the datasets. 

In Table [11] we present the constraints on Yhc, rUe, 
and the rest of the cosmological parameters for all of our 
datasets. For the CMB dataset, Yho = 0.296 ± 0.029, which 
is almost 2 a higher than the standard value. The constraint 
on me is me/(me)o 
correlation function 



0.9Mt°o°l°8- By adding the CMASS 



factor of three (me/(me)o = 0.985t!5:^2? 

) but the bound 



we tighten the bounds on irie by a 
)o = 0.985tH^?) 
on Yue is unchanged. When all datasets are considered in 
the analysis, the constraints are Ihc ~ 0.302 ± 0.029 and 
me/{me)o = 1-0061q oi4. The Ho prior shifts Yuc towards 
slightly higher values, which are more than 2 a above the 
standard value of 0.24. On the other hand, rrie is consistent 
with no variation at 1 a, for all of the datasets. 

Regarding the values of the remaining cosmological pa- 
rameters, we find that there is a slight tension at the 1 a level 
for the parameter ris, which is higher than its ACDM value 
for the CMB+CMASS dataset, and for the parameter erg, 
which is higher than its ACDM value for the full dataset. 
In the rest of the cases, the parameters are all consistent 
within 1 a with their ACDM values. 



Table 10. The marginalized 68% allowed regions on the cos- 
mological parameters of the ACDM model, adding the variation 
of the fine structure constant a and Yfje as a free parameter, 
obtained using different combinations of the datasets. 
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Figure 16. The marginalized posterior distribution in the me 
- Yho plane for the ACDM parameter set extended to include 
the variation of and a variable primordial helium fraction 
Yfjc- The dot-dashed lines show the 68% and 95% contours ob- 
tained using CMB information alone. The dashed lines corre- 
spond to the results obtained from the combination of CMB 
data plus the shape of the CMASS ^{s). The solid lines in- 
dicate the results obtained from the full dataset combination 
(CMB-l-CMASS+BAO+Ho)- 
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Table 11. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, adding the variation of 
the electron mass me and Yhc as a free parameter, obtained using 
different combinations of the datasets. 
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4 CONCLUSIONS 

In this paper we have presented new constraints on the vari- 
ation of the fine structure constant and on the electron mass 
using the latest CMB observations, and the full shape of the 
(spherically averaged) redshift-space correlation function of 
the CMASS sample of galaxies, drawn from the Data Re- 
lease 9 (DR9) of the Baryonic Oscillations Spectroscopic 
Survey (BOSS). Recent BAO and Ho measurements were 
also considered. We have studied the degeneracies between 
these constants and other cosmological parameters, such as 
the dark energy equation of state, the neutrino mass, the 
effective number of relativistic species, and the primordial 
helium abundance. The main results can be summarized as 
follows: 

(i) In the case of variation of only a, our bound is 
a/ao = O.9957I0 0042 5 consistent with no variation at the 
2 a level (almost at 1 a). The constraints from CMB data 
alone are slightly improved when CMASS and other cosmo- 
logical datasets are included. 

(ii) When only me is allowed to vary, the bounds are 
highly improved when additional datasets are considered 
in the analysis. Our best estimate including all datasets is 
me/(me)o = 1.006^0 013 5 consistent with no variation of nie 
within 1 a. The CMASS dataset improves largely the CMB- 
only constraints. 

(iii) When both a and me are allowed to vary, the con- 
straints on a do not improve with the addition of new 
datasets, while they do for me. With each dataset addi- 
tion, the variation of the constants becomes more corre- 
lated. Our final bounds, using the complete dataset, are 
a/ao = Q.mQlt^omfi and me/(me)o = 1.028 ± 0.019. Both 
limits are consistent with no variation of the fundamental 
constants within 2 cr. 



(iv) When we study the variation of a taking the value 
of wde as a free parameter, the allowed region for a hardly 
depends on the addition of further datasets. The constraints 
on wdb, however, are noticeably improved. The bounds for 
the complete datasets are a/ao = 0.9915 ± 0.0048 and 
wdb = — 1.20±0.13, which deviates from the standard value 
jude = — 1 at 1 (t. 

(v) When both me and u;db are allowed to vary, the 
CMB dataset cannot place tight constraints on both quan- 
tities. The inclusion of CMASS correlation function helps 
to improve the bounds. Our results for the full dataset are 
me/(me)o = 0.996 ± 0.029 and wdb = -1.12 ± 0.23. 

(vi) When we allow for a non-negligible contribution of 
massive neutrinos to the dark matter component, our result 
is a/ao = 0.9978lo.oo45 - and < 0.043 at 95% CL for 
the full dataset. The constraint on the sum of the neutrino 
masses is rrii, < 0.46 eV at 95% CL. 

(vii) In the case of joint variation of me and the 
bounds for the full dataset are me/(me)o ~ 1.035 ± 0.021, 
and < 0.086 at 95% CL. The degeneracies between these 
quantities and Ho shift the constraints to higher values when 
the HST Ho prior is included. 

(viii) When A'cS and a are allowed to vary, the values of 
A^ofi are consistent within 1 a with its standard value for 
the CMB and CMB-^CMASS datasets. When the HST Ho 
prior is added, the value of A'^off is much higher, being slightly 
consistent at 2 a with its standard value. The constraint on a 
is a/ao = 0.9967± 0.0042, fully consistent with no variation 
at 1 a. 

(ix) When A'ofi and rrie are allowed to vary, the value of 
Acff is more than 1 a higher than its standard value, and 
the value of me is consistent with no variation at 1 a, for all 
of the datasets. 

(x) In the case of joint variation of Vho and a, there is 
a strong degeneracy between these quantities, for all of the 
datasets. For the CMB and CMB-f CMASS datasets, Yhc is 
consistent with the standard value of 0.24. When the Ho 
prior is added, Yhg is almost 2 a higher than the standard 
value. For all of the datasets, there is no variation of a at 
the 1 a level. 

(xi) In the case of joint variation of Ihc and me, there is 
no correlation between these quantities for any dataset. The 
bounds on Yue are similar for all of the datasets, a value 
around 2 a higher than the standard one. For all of the 
datasets, me is consistent with no variation at 1 ct. 

(xii) We present new bounds on A'^eff and in the 
case of no variation of fundamental constants. For the 
full dataset, in the ACDM-|-Acff model, the constraint is 
Ncs = 3.86lo4o. In the ACDM-|->He model, the constraint 
isFHc = 0.30li°;°29- 



The analysis carried out in this paper is based on the 
first spectroscopic data release of BOSS. Future data releases 
will provide even more accurate views on the LSS clustering 
pattern by probing a larger volume of the universe. Together 
with Planck satellite data to be released in early 2013, this 
will enable to put more stringent constraints on the variation 
of fundamental constants. 
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Table Al. The marginalized 68% allowed regions on the cosmo- 
logical parameters of the ACDM model, obtained using different 
combinations of the datasets. 
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APPENDIX A: CONSTRAINTS ON THE ACDM 
MODEL 

In Table lAll we present our constraints on the cosmological 
parameters of the ACDM model, for better comparison with 
the results obtained in t his paper. These are consistent with 
the constraints given in ISanchez et aP (|2012l ). 
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